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ABSTRACT 


In the present analysis approximate natural frequencies 
of rotating orthotropic discs have been obtained. Ritz 
method is used assuming approximate mode shapes. Radial 
variation of thickness and density are also considered. 

Effect of clamping radius, variation of thickness, variation 
of density and variation of rigidity parameter (ratio of the 
Young's modulus in tangential direction to the Young's modulus 
in radial direction) are studied. Special cases of rotating 
isotropic discs and stationary orthotropic discs have been 
compared with previous work and found that the frequencies 
obtained from the present analysis are well in agreement 
with them. 

It is found that frequency of rotating orthotropic 
discs increases with increase in clamping radius. As the 
disc becomes thicker at the hub or heavier towards the 
periphery frequency increases. The effect of increasing 
the rigidity parameter is to increase the natural frequency* 



CHAPTER 1 


INTRODUCTION 

1.1 General Introduction: 

Vibration characteristics of circular spinning discs 
are of interest in design, since circular discs are commonly 
used structural elements as in turbines, flywheels, saw 
blades, optical and radar reflectors for space vehicles, solar 
sails etc. It has been observed that while a turbine disc 
or a flywheel is kept spinning it exhibits transverse vibra- 
tion. Research work so far made in this connection, have 
shown that such vibrations of disc are largely responsible 
for the disc failure, fatigue of the metal and gradual deve- 
lopment of cracks on the disc or wheel. Such vibrations at 
certain speeds, become veiy pronounced and thereby cause 
considerable additional bending stress which may result in 
fatigue of the metal and gradual development of cracks. 
Sometimes vibrating discs majr even come in contact with 
adjacent parts and lead to the failure of the whole structure. 
The necessity for the design of efficient, lightweight 
structures has led to the use of special alloys and composite 
materials which may have directional properties. Similarly, 
design may demand for a particular type of use or to improve 
mechanical strength and properties of circular plates in 
desired directions, incorporating stiffness radially or 
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circumf erencially or in both* Stiffened plates oan be 
idealized as orthotropic plates. Thus a study of the 
vibration characteristics of a rotating isotropic and ortho- 
tropic plates becomes essential. 

■~"1.2 A. Brief Survey of literature: 

\ , 

1,2.1 Isotropic Discs. 

Early investigators found that the solution of the 
general case posed considerable mathematical, complexity and 
elected to study separately the cases in which the rotational 
effects were negligible and the cases in which bending 
stiffness effects were overshadowed by rotational effects. 

The case in which rotational effects are of primary importance 
led directly to the study of spinning membrane discs, lamb 
and Southwell [1] discuss the problem of vibration of a 
spinning circular membrane - a disc with negligible bending 
stiffness - and it is shown that the mode shapes are Jacobi 
polynomials, Southwell [2] extends the analysis to the study 
of the transverse vibrations of a spinning membrane clamped 
to a rigid hub in such a manner that radial but not transverse 
displacement may occur. The solution to this problem 
necessitates finding of hypergeometric functions. Simmonds[3j 
considers the same problem and points out that Southwell’s 
paper contains an error since it overlooks the facts that the 
hypergeometric equation cannot have two independent solutions 
invoicing only teras containing the hypergeometric series 

i Y»x) , if y is allowed to take on negative integral value 
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Bulkeley and Savage [4] have studied the centrally 
clamped membrane for the case of axisymmetric vibrations. 

Their results are for various degrees of partial clamping, 
including the fully clamped configuration. Transverse 
vibrations' of centrally clamped membrane with Coulomb 
friction at the disc collar interface, are shown by a scale 
factor, to be equivalent to those of a membrane fully built 
into a hub, or clamped between frictionl ess collars. 

Simmonds [5J has studied the fully clamped case for axisym- 
nctric vibrations, in which equation of motion reduces to 
the hypergeometric diff erential equation. 

Mote [ 6 'J* has on ployed a Bayleigh-Hitz technique to 
study the approximate vibration characteristics of centrally 
clamped variable thickness discs. He investigated the 
natural frequencies of transverse vibrations taking into 
consideration rotational and thermal in-plane stresses as 
well as purposely induced initial stresses. He found that 
initial stresses can significantly raise the minimum disc 
natural frequency throughout a prescribed rotational .and 
thermal environment. The fundamental node of disc vibration 
is one of zero nodal circles and either zero, one or two 
nodal diameters, depending upon the disc geometry and the 
rotational- thermal environment. Mote further extended his 
work [ 7 J to incorporate the transient temperature distribution, 
and associated thermal stresses. An analytical rul 0-0 f- thumb 
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involving the calculation and. test of a non-din ensional 
number is proposed for predicting conditions when the 
transient regime nay be significantly more severe than the 
steady state. He has also investigated [3j the possibility 
of increasing the fundamental frequency by purposely inducing 
thermal membrane stresses in centrally clamped, poripherilly 
free rotating annuli. Practical application of both these 
methods [6,8j is limited to the ceases where minimization of 
disc thickness is of paramount importance. Optimization of 
the fundamental frequency and limitations of the above 
methods are also discussed [8,9]. 

Eversman [I0,llj has given the exact method for 
determination of natural frequencies of transverse vibrations 
of a clamped spinning membrane disc. La,ter on, Eversman and 
Dodson [12] treated the problem of transverse vibrations of 
a spinning centrally clamped circular disc, for which Mote [6j 
has given an approximate solution, Their approach is exact 
to within the accuracy of the numerical computations, requires 
no appropriate choices of the nodal functions, and is free of 
usual uncertainities of approximate methods, Barasch and 
Ghen [13] have also investigated the above case of rotating 
disc, the equation of motion is solved by reducing the fourth 
order equations of notion to a set of four first order 
equations subject to arbitrary initial conditions. A modified 
Adams’ method is used to numerically integrate the system of 
differential equations. 
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1.2*2 Experimental Analysis. 

Not much experimental analysis on the rotating discs 
is reported in the literature. Campbell [14 J investigated 
experimentally various forms of vibrations and waves which 
may exist in steam turbine disc wheels and proposed 
procedures necessary for the protection of steam turbine 
bucket wheels from transverse vibrations, krauter and 
Bulkeley [l5j have experimentally observed the free trans- 
verse vibration of a centrally clamped spinning membrane 
discs and compared axisynnetric and asymmetric frequencies 
with theoretical predictions [ 3 ]. 

1.2.3 Ortho tropic Disco. 

Recently use of composite materials and stiff fining 
of circular plates has led to the investigation of the 
characteristics of orthotropic circular plates. Circular 
plates with, different elastic properties in radial and 
circumferential directions have been called by various 
names as cylindrically aeolotropic, orthntropic, polar 
orthotropic or anisotropic by different investigators. 
Minkarah and Hoppnann [l6j consider flexural vibration of 
cylindrically aeolotropic circular plates sis cun ing 
symmetric vibrations. The theory is exemplified by an 
analysis of the vibration of a circular plate with circular 
stiffness for which the effective elastic compliances had 
been determined in previous investigations. 


They also 
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discuss the approximation of equivalent homogeneous unifomly 
thick aeolotropic plates for stiffened plates. Dzialo and 
Hoppaann [17 J extended the above problem to consider the 
effect of rotaiy inertia. They found that the frequencies 
are higher than the actual if rotary inertia is neglected. 
Pandalai and Patel [l8j have also considered the above 
problem and determined natural frequencies of orthotropic 
circular plates. 

Ghosh [l9j in vestige, ted vibration of a rotating 

m 

aeolotropic elastic circular disc of uniform thickness. He 
has given a mathematical analysis fir determining frequencies 
and nodes of vibration. He further extends his investigation 
[20j to take into consideration disc of variable thickness. 

In his analysis, he neglected the effect of bending stiffness, 

Hanaiah and Vijaykunar [21,22j have used Hayleigh-Bitz 
method with simple polynomials as admissible functions for 
obtaining natural frequencies of transversely vibrating polar 
orthotropic annular plates. Estimates of natural frequencies 
corresponding to nodes with one as well as two nodal diameters 
are obtained for the nine combinations of clamped, simply 
supported and free edge conditions and for various values of 
rigidity ratio and hole sizes. Based on the variation of 
eigenvalue parameter with rigidity ratio, the frequencies of 
these nodes as well as tho axi symmetric no^es have been 
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expressed by means of simple formulae in terms of rigidity 
ratio and corresponding mo^es in the isotropic case. 

1.3 Present Work : 

The present work deals with the determinations of 
natural frequencies of orthotropic rotating discs. The effect 
of clamping radius, thickness variation and the density 
variation are determined. The effect of rigidity parameter 
( The ratio of Young’s modules in Q direction to the 
Young's modulus in r direction) on the Natural frequency is 
also ahalysed. Approximate natural frequencies are obtained 
by using Ritz-method . Specialising the above results to the 
case of an isotropic disc, comparison is made to existing work. 



CHAPTER 2 


PROBLEM FORMULATION AND SOLUTION 

The necessary formulation and solution for the 
calculation of natural frequencies of an orthotropic circular 
disc of variable thickness and density, rota.ting at a constant 
angular velocity is presented in this section. Expressions 
for the potential energy in bending and potential energy due 
to rotation (that is, due to the inplane radial and 
circumferential stresses) and the kinetic energy of the disc 
are given. In-plane stresses for the above orthotropic, 
rotating disc is presented. Assuming that the disc is under- 
going harmonic oscillations, suitable mode shapes have been 
assumed and expressions for maximum potential energy and 
kinetic energy are obtained. Hits method is used to obtain 
the natural frequencies. Frequencies for different plate 
geometries, variation of density and thickness and orthotropy 
in two directions are obtained. 

2.1 Plate Geometry: 

Consider a thin orthotropic elastic disc of outer 
radius ‘b*, clamped at radius ’a f , rotating uniformly with 
angular velocity fl about the normal axis passing through the 
centre and perpendicular to the centre plane of symmetry as 
shown in Fig. 1. 
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Variation in the thickness of the disc has been 
chosen in the form 


h 



( 2 . 1 ) 


where, h = thickness at radius r 

h Q = thickness at periphery 

b = outer radius 

|3 = thickness parameter. 

Similarly variation in the density of the disc is 
assumed to be of the form of 


„ m 

o = P 0 (f) (2.2) 

where, p = density at radius r 
P Q = density at periphery 
m = density parameter. 


2,2 Boundary Conditions; 

The geometric dispL cement boundary conditions at 
the clamping radius r = a a.re 

(i) Transverse displacement * w’ is zero, i. e. 

w = o at r = a (2.3) 

(ii) Slope is zero, i.e. 


= o at r = a 


(2.4) 
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The following boundary conditions are used to determine 

the radial and tangential stresses due to rotation. 

(iii) Badial displacement ’u* is zero, i.e, 

u = o at r = a (2.5) 

and 

(iv) Badial stress is zero at the periphery, i.e. 

or = o at r = b (2.6) 

r 

2.3 Assumptions i 

This analysis is based on the following assumptions: 

(i) Conditions of plane stress. 

(ii) The thickness and density distribution is symmetrical 
with respect to both the axis of the disc and the 
midplane of the disc. 

(iii) The influence of shear stress and rotary inertia is 
negligible. 

(iv) The disc material is polarly orthotropic. 

(v) The stress-strain relations follow a generalised 
Hooke’s law with the principal axis of orthotropy 
coinciding with the principal axis of the disc. 

2.4 Potential Energy due to Bending: 

The potential energy for the plate in bending can be 

written as [22j, 
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v B . * r w*? a ♦ 2 ^(4? - i 


o a r Sr* 


Dr 


r 3 r 2 ar r ae 2> 


+ + 1 -jiw ) 2 + ? D 1_(1 JL-S) 2 lrdrd6 

„ 2 V 3r r ao t ' + ^ r6 3r l r 3 «' JrclrcW 

(2.7) 


36 


where 


E h 3 
X* 

D = * % t flexural rigidity in radial direction 


Dg - 


12(l-v r v Q ) 


V 3 

12(1- v r v g ) 


( 2 . 8 ) 

t flexural rigidity in tangential 


direction 


(2.9) 



v r D 6 ”^6 D r 


( 2 . 10 ) 


D r© 55 Gh.Vl2 t , shear rigidity 


(2.11) 


f(B r B 0 ) 

2[l-/( MpVg) j 


(2.12> 


E^ » Young’s Modulus in radial direction 
Eg = Young’s Modulus in tangential direction 
v y * strain in tangential direction due to the unit strain 
in radial direction 

Vg = strain in radial direction due to unit strain in 
tangential direction. 

2.5 Potential Energy due to Rotations 

Potential energy due to rotation can be written 


as [23j» 
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2 IOC *fc) ^ pljA a a JE o 

v a = 4 J 1 [to r C^j) 2 + te e ( 3 B |) 2 Jrdrde (2.13) 

O Q. **■ 

where, cr_ — radial stress 
r 

and = circumferential or tangential stress. 

2.6 Kinetic Energy i 

Kinetic energy of the disc can be written as [23 J, 

2% b 2 

1 = i j J P h(-^-~ £) rdrd© (2,14) 

o a '' 

2.7 Expressions for in-plane Stresses; 

The expressions for the in-plane stresses of a 
rotating orth^tropic disc with variable thickness and variable 
density has been obtained by Reddy and Shrinath [24]. The 
method is giv en here b ri ef ly f 

For a spinning disc rotating at an angular velocity ft, 
equation of equilibrium to be satisfied by the stresses is 


~(hrcr r ) - hff£ + hP^ 2 r 2 = 0 

The strain disp.Tac ement relations are 

_ du _ u 

e r * dr r ~ r 

and hence the strain compatibility equation is 


(2.15) 


( 2 . 16 ) 


v dr ve r e 0 ; 


0 


(2.17) 
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The stress-strain relations for the polarly 
material are 


_ fr 

e r ~ S r - E 0 


'© “ E, 


£ 

b 


^r „ 

E a r 
r 


With Maxwell’s relation 


v, 


G 


E, 


G 


E 


ortho tropic 


(2.18) 

(2.19) 


( 2 . 20 ) 


The equation of equilibrium can be satisfied by a function P 
such that 

P = h o p and h cr Q = + hp ft 2 r 2 (2.2l) 

Using Eqs.(2.18 _ 2.20) ajLongwith the relations (2.1) and 

( 2 . 2 ), 


h » h 0 (r/b) - ^ 
and p as p 0 (r/b) m 

in Eq.(2.17) we get, 

43E + a? _’fk , _ K h . o! m-p+l h m-p 
*^ras:n F = S h o p o a r b 

where, - 1+P 
^ = X2 + 

K^ ~ — ( 3+JO+ ^q) 


( 2 . 22 ) 
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Using Eqs.(2.15 - 2.22) radial and circumferential stresses 
been 

have/ obtained as discussed by Reddy and Shrinath [24,3* ^ or 
the disc clamped rigidly at the hub, the boundary conditions 
are; 

(i) displacement is zero at clamping radius, i.e. 

u = o at r = a 

and 

(ii) radial stress is zero at the periphery, i.e. 

a =o at r = b 
r_ 


Eor the above conditions stresses are given by 


P an+P-1 p_ a 0 +(3-l 0 

= p„ C 3 [|l(r/a) 1 - ^(r/a) 2 + (r/b) 24ffl J 


(a. 23) 


a, p. 


a Q - 




tt 1+ P-l a 2 |3 2 a 2 +p-l 


3 P, 


P; 


2 +m 2 +m 

+ ( 3-Hn-P) ( r/b) } + (r/b) J 


(2.24) 


where, 

a i = - + V(^ P 2 + x 2 + ^ p) 
a 2 = - -|rP - VX4 P 2 + X 2 + Vq p) 
X 2 — 

-( 3 + m + v Q ) 

Q = ~. TW imi - ,. LJ v y 

^ (9-X 2 )~P(3+m+v c ) + m(m+6) 


(2.25) 

( 2 . 26 ) 

(2.27) 

(2.28) 
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“X 2 3+m -a 0 

) (a/b) + (a 2 - v 0 )(a/b)^ ( 2 . 29 ) 

v Q - X 2 3+m-a, n 

P 2 = ( 3 ^+^—) (a/b) 1 + (a x - v^Ma/b ) 13 ( 2 . 30 ) 

and 

P 5 = (a x - v^) (a/b) ^ - (a 2 - v Q )(a/b)' ^ (2.31) 

It may be noted that the stress and displacement 
expressions are not valid for combinations of X 2 , v^, m 
and p which make 

(9 - X 2 ) _ p (3 + m + v^) + m(m+ 6 ) =0 (~.32) 

2.8 Assumed Mode Shapes i 

Natural frequencies of the disc have been obtained 
by the Sitz method which is b sed on the principal of 
minimum potential and involves determination of maximum 
potential energy and maximum kinetic energy [2 5 j» as well 
as assumption of appropriate mode shape function, which 
satisfies atleast the geometric boundary conditions and 
approximates the actual modes of vibrations. 

Assuming that the plate is undergoing harmonic 
oscillations, the lateral deflection *w T can be written as 

w — W^(r, ^) cos ^ 'fc (2.33) 

and 

¥^(r,©) is separable as 
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W 1 (r, 0) = W(r) cos n © (2.34) 

where n is the number of nodal diameters. Note that the 
geometric boundary conditions (2.3) and (?.4) depends only 
upon the function, W( r) . The mode shape W(r) in the radial 
direction is chosen as [6j 

VJ(r) = a Q (r-a)^ + a-^(r-a)^ + a2(r-a)^ + a^(r~a)^ + a^(r-a)^ 

+ a^(r-a)^ (2.35) 

where a Q , a-^, a^, and a^ are unknown constants. 

2.9 Maximum Energies: 

Substituting for ' w f from El. (2.33) in the expressions 
for potential energies, Els. (2.7) and (2.13), we get 


2jc b 3*W- p 231 3 2 3L 3W_ _ 3%_ 

V_ = 4-cos 2 «» t } / [D ( -) + — — — ( + - -) 

B ‘ n O a r r 3r r 3e* 


'1,2 . 2 \ 


D„ * W. , 9* V, 2 , , 3 «, ,2 

4 . at 4 . i — ) 4 . pin f - S-- (1 — ui) } lr*drd© 

r a' Sr r ; “ J r© 1 3 r'r 3 Q } jrara 


(2.36) 


2tc b 3 \L 2 3 W, 2 

= 4. cos 2 c^t J J { ho r ( -g“) + ha^-^jpg) } rdrd© 

o a 

{ 2 , 31 ) 


The potential energy of the plate is maximum when 
the displacement is maximum, which occurs when cos w^t equals 
unity. This is true when 
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t = XLjTi (el^s o, 1, 2, 3, ...} (2.38) 

Using these values of a> n t, the maximum potential energy 
due to bending and rotation are obtained as 


B max 


= * 


2n b 3 2 17, 9 2D, 32 % 9 W, , 3 2 1 L 

; / [dj i) 2 + i -4) 


r 9 
o a 9 r 2 


a r 


2 3 r r 


b q . 3 \ 1 

+ ~(-r™ + - 

r 2 3 r r 


d z 17 , o , , 3 V- o 

i -jgi) } 2 Jrdrae 


3 4 


(2.39) 


B. max 


271 b 3 W, ~ , 3 17, 0 

i S I {h + ho G (i Yq) > rdrdG 

00 ~ " 


(2.40) 


Similarly, from Els. (2.33) and ( 2 . 14 ) kinetic energy can be 
written as 

271 b 2 

T = 1 sin 2 w t « * | j p h W, rdrdQ (2.41) 

n n 0 b - 1 


The kinetic energy is at a maximum, when the velocity of tho 
plate is maximum, which occurs when w(r,4,t) is zero. This 
will be true if cos w n t = 0; thus 

w n t = (a x + -3r)it ( 2 ^ = 1,2, 3,4, 5,...) (2.42) 


Substituting these values into (?.4l) 
for maximum kinetic energy as, 


w 


2 it b 

- +1 J ****** 


a 


we obtain the expression 


T 

max 


(2.43) 
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Substituting Eq.(2.34) in Sis. (2.39), ( 2 . 40 ) uni (2.43), we 
get 


ma * - * 7 J too-- <v-0> 2+ -£*> 


D , 


+ -S(4? - ?! V )2 }+ sin 2 n© ' 3 * * 2 


n 2 


,2 3 r r 


i-fi - p‘]rdrde 


(2.44) 


V 


2% b 


R max 


\ J J ( ha^cos 2 n©(*y~) 2 + hcr a n 2 sin 2 n© ( 7/ r) 2 }rdrd© 


o a 


©‘ 


(2.45) 


and 

2tc b 

T max = ^ w n i" f p hcos 2 n©7 2 rdrd© (2.46) 

' o a 


Therefore, since 

2% 2tl 

j sin 2 n0d0 = J cos 2 n©d© = % (2.47) 

o o 


we get 


v = | j b [D (43 )2 + ! 2 k f !( 5 S 

B max 2 j L r . « " 


a 


dr 2 


P 2 

. - ~ w) 

flr a dr r 


+ ^ m . Sl W )2 + i2s«H(42 - 2)2j rr , r ( 2 .48) 

„* dr r „ 2 dr r J ^ ' 


V . 


B. max 


% 

2 


a 


Chc r (f) 2 


+ hcgii 2 ( W/ r ) 2 J : 


and 

_ It 2 

T =fu 
max 2 n 


b 

J p h¥ 2 rdr 
a 


(2.49) 


(2.50) 
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2.10 Uon-dimensionalization r 


To express energy expressions (2.48 - 2,50) in 
non-dim ensionalised form, we set 


r =s hr 

( 2 . 51 ) 

Substituting Eqs.(2.l) and 

(2.2) for h and P , and 

expressing D r , D Q , and 

IL^ in the following form. 

D r " \ 

(2.52) 

** (3 

= v 

(2,53) 

_ - --3P 

D r© * r 

(2,54) 

- _-3P 
\ * *L r 

(2.55) 

where, 

= 

(2.56) 

12(1- v r \> Q ) 

5 !k.-8 

w 12(1- v r v Q ) 

Gr h 3 

D r0 ^ 12 

(2.57) 

(2.58) 

\ = \ \ - v e 5 r 

(2.59) 


The expressions for maximum energies can now he written as. 
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\ pr 2 \rd*U n* t 


v. „ -SL / rii rsjjjh* + —ir: 

B max 2 v2 /, ‘-A r \ 'gJ + - '•.-2 _ 

b z a/b dr r dr r 


'/) 


- 3<« ~ « 2 - ^ dr 


- r dr r 


r dr r 


V R max * 2 Ai o — r 


and 


(2.60) 


? /, [cr^(^) 2 + n 2 a G (V/r) 2 ] r ^ dr (2.6l) 

dr 


_ o 1 - 1+m-p 

T = r b f d h W* r ^r 

max 2 p o o r " r 


(2.62) 


Setting 


r - a as x. 


(2.65) 


where a = a/b 


(2.64) 


Now Bq.s.(2,60) # (2.6l) and (2.62) can be written as 


% \ 1 7 a r ad 2 ¥ x2 , \ 2 

'B max - 2 b2 J U dx 2 } 


^(® - -a!— ¥) 


5^ (x+ a) dx a c1jc ( x+a) 


'© /d¥ _ n 2 ,,\2 , 4D r© n 2 


D r (x+a) 


(J 

-\2 '‘dx 


( x+a) 


■V/ ) ^ + 


( x+a) 2 


( 


dV 


dx (x+a) 


w a 1-3P 

j ( x+a ) dx 


(2.65) 


Similarly, 


7 h Bax - K °r ( S> 2 + °e 

0 A . (2.66) 
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-Ed. 


E 


1-a 


max 


f b 2 p Q h Q J 1 ./ 2 (x+a) 1+m_ J dx 


(2.67) 


Qe assumed mode shape (2.35) takes the form 


¥( r) = a x 2 + a, x 3 + a,x^ + a^x^ + a.x^ + (2.68) 


He the same manner expressions for radial and tangential 
stresses can be written as 


cu +P-1 
h (x+a) 1 

‘ ’o -5 ^ 0 - 

3 (a) 1 


^2 


a 0 +f3-l 
(x+a) 2 


3 % ( S ) “ a *?- 1 


+ o 3 (^) 2+m J 


(2.69) 


a_Ed 


8*a = 


_ a, +(3-1 
a 1 p 1 (x+a) 

p ft b 2 TCL — afe— air 

0 Lh 3 P, a,+p-l 

3 (5) 1 


a. 


__ CE2+P-I 
^2^2 ( x+ ^) 

3 a 0 +|3— 1 

3 (5) 2 


+ Cj( 3+m+P) ( 3c+-a) 2+m + C;i+a) i " r!a ] 


, — \ 2+m • 


(2.70) 


f*«Llowing non-dimens ionalized terms are defined 


p -^T * — t non-dimensional radial stress (2. 71) 




p 0 0 b 2 


, non-dimensional tangential stress (2-72) 
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XI' 2 p Q h 0 b 4 

£2 = V" ( — — — — )» non-dimensional speed (2. 73) 

8 Sr 

— P n h n b 4 

= V" ( — - — - ), non-dimensional frequency (2*74) 

5 

r 

2*11 Determination of Natural Frequencies : 

Substituting the expressions for the mode shape 
from Eqn. (2*68) and for stresses from the Eqns*(2.69) and 
(2*70) in the expressions for maximum energy Eqns .( 2*65-2* 67 ) 
we obtain 

W ir ^ S r f i (a i- B e • V v e < -■»•“) ( 2 -? 5 > 


Bmax 2 0 


P Q 8 2 b 2 ip^i’S©’ V© »a,r,p,m,n) (2*76) 


'"2“ ' &q Pq u n S » e ? r j p > m ) 


(2*77) 


Applying Ritz— method for obtaining ■frequencies » we may write 


aUL =0 (i = 0 , 1,2 5 ) 


(2*78) 


where 


tr > *tr fp 

Bmax Rmax max 


(2*79) 


from Eqns. (2*75-2*77) we obtain. 


9 7 Bmax h Q p 0 ** b 3 Y Rmex 

3 5 3a i 


h 0 p o“n 1,4 


3 a i 


(2*80) 
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where [V^], [V R 3 and [T max 3 are 6x6 matrices 
and 



Eqn. (2.81) can also be written as 


( 2. 81 ) 


(2.82) 


[K] (X } = «2 [m] { X> (2.83) 

where 

[K3 * [ T b 3 + ST J [V r 3 (2.84) 

and 

W = [ Vr 1 <2- 85) 

-1 

Bqn. (2.83) can be premulti plied by [M] to give 

[M]" 1 [K] { X } = m ^ <X> (2-86) 

and can be reduced into an eigen value equation, 

[[A]- m ^ [I) ) { x } =0 (2.87) 

Natural frequencies are obtained from the roots of 


e qn . (2*87). 



CHAPTER 3 

RESULTS AND DISCUSSIONS 

Natural frequencies of orthotropic discs have been 
calculated numerically from Eqn. (2.87). The effect of variation 
of clamping radius ^ , thickness parameter {3, density parameter 
m, rotational speed Q , rigidity parameter X^ and poisson's 
ratio have been studied ^Frequencies have been obtained forO to 
5 nodal circles and for 0 to 4 nodal diameters. Frequencies 
plotted in Figs. 2 » 7-8, 10-14 and in 16-18 are for zero nodal 
circle and zero nodal diameter mode. The parameter D^q/D^, was 
fixed at 0.35 and Poisson's ratio Vg was fixed at 0.3 
(except in the case, in which the effect of variation of 
is studied). Computations were carried out with double precision 
arithematic for all type of discs. Comparison have been made 
with available results for isotropic and orthotropic discs. The 
frequencies obtained in the present analysis is in good 
agreement with previous work. 

3.1 Effect of Variation of Clamping Radius 

p 

Natural frequencies of a stationary, isotropic (A =1.0) 
disc of constant thickness and constant density have been 
obtained, for various clamping radius ^ and compared with 
Mote [6], (Table -1). It is found that present frequency 
estimation compares well with that of Mote. Frequencies obtained 
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for above case have also been compared, with exact natural 
frequencies given by Southwell [23. Approximate and exact 
natural frequencies for various clamping geometries and at 
different nodal circles and diameters are shown in Table 2. 

For ^ < 0.7 the present approximate values of the natural 
frequency compares well with the exact values of southwell 
(Fig. 2). For | > 0.7 the error is large , however, this range 
is not of practi cal interest . Fig. 2 shows plot of clamping 
radius vs. natural frequency at zero nodal circle, zero nodal 
diameter for present, Mote's and Southwell's case. 

The variation of natural frequencies with clamping radius 
^ is shown in(Fig. 3) in the case of a disc of constant 
thickness and density, rotating at a speed of sT = 20. The value 

p 

of the rigidity parameter A is taken to be 0.25. Figs. 4 and 

p 

5 are plotted for the case when X =1.0 and 4.0. Frequencies are 
plotted for various nodal circles and nodal diameters. From the 
above analysis it is found that the natural frequency increases 
as the clamping radius increases from .1 te .8. 

3.2 Effect of Variation of Thickness Parameter : 

Natural frequencies of a stationary isotropic disc 
( X 2 = 1 . 0 )kave teen obtained for |3=0. 9 and are compared with the 
results of the Mote [63 (Table 3.). The comparison is very good. 
From table 1 and 3 it can be seen that frequency is high for 
{3=0. 9 than for p=0.0 (constant thickness). Turbine discs can be 
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dealised as a disc whose thickness decreases hyperbolically 
n radial direction (Pig. 1-b) . it is found that as thickness 
arameter p increases from 0 (constant thickness disc) to a 
ositive value (that is the disc is becoming thicker at the 
ub), frequency also increases as can be seen from Pigs. 7,10,11 
nd 15* On the other hand if j3 increases negatively (that is the 
isc is beo-oming thicker at the periphery) the frequency 
3 creases. Pig. 16 shows variation of frequency with the nodal 
Lameters for different values of p. 

■ 3 Effect of the Variation of Density Parameter ' 

Use of composites may cause radial variation of density, 
ieping this in mind natural frequencies of disc having 

p 

= 1 and 4.00 respectively with different density 

rameter m, have been determined and plotted in Pig. 8. As 
e value of parameter m increases ( that is density is more 
wards the periphery) frequencies decrease for X=' 

00 but for X^= 4.00 it slightly increase. 

4 Effect of Rotational Speed: 

Effect of rotational speed on the natural frequencies 
isotropic and orthotropic discs have been studied for various 
sc geometries and density parameter m. Frequencies of an 
^tropic rotating disc of constant thickness and density have 
m compared with the work of Baras ch and Chen [133 and 
ithwell [2]* Frequencies are plotted for 0, 1 and 2 nodal 
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circles and for zero nodal diameter (Fig. 9), which agree well 
with the frequencies obtained by -Barasch and Chen and is slightly 
greater than the Southwell’s values. Figures 10 to 14 show the 
variation of frequency with the rotational speed with X 2 as a 
parameters for various values of {3 and m. is the rotational speed 
increases frequency also increases. This is because of in-plane 
stresses due to rotation. Increase in speed increases in-plane 
stresses which inturn increases stiffness and hence the frequency. 
Fig. 16 shows the variation of frequency of a membrane with rotation. 
3*5 Effect of Variation of Rigidity Parameter: 

Frequencies of stationary orthotropic disc of constant 
thickness and constant density for different clamping geometries 
have been obtained and compared with the frequencies obtained 
from the approximate formula given by Ramaiah and Vi jay Kumar 
[22) in Table 4. Frequencies obtained from the present analysis 
are in good agreement with those obtained from approximate formula. 

Figs. 10 to 14 shows variation of frequencies in the case 
of a rotating disc for different values ofX 2 . It is seen that 
at lower speed the frequency increases as x 2 increases but at 
higher speeds the frequency decreases with increasing x 2 • When 
a rigid rotating disc is executing bending vibrations, it is 
controlled both by bending stiffness and tension due to 
in-plane stresses. Since the tensions in the disc are proportional 
to the square of the angular velocity, it is clear that at low 
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speed the main controlling factor is the bending rigidity. 

Whereas at high speeds the controlling factor is mainly the 

in-plane stresses. The effect of increasing X^ is to increase 

the natural frequency. The stresses due to rotation decreases 

with increase in the rigidity parameter X^ (which is shown 

in Fig. 6) and consequently the natural frequency decreases with 
2 

increase m x • This explains why the natural frequency 
increases with increasing X at lower speeds and decreases 
with increasing at higher speeds. In the case of membrane 
in which the bending rigidity is zero, the variation of natural 
frequency with speed at different X^ is shown in Fig. (17)* 

3.6 Effect of Poisson’s ratio ; 

Fig. 18 shows the variation of frequency with Poisson’s 
ratio Vq for a constant thickness , constant density isotropic 
( X^=l) and orthotropic ( X^=4 ) discs rotating at speed ^=15* 



CONCLUSION 


In this investigation the effect of clamping radius , 
thickness parameter, density parameter, rotational speed, 
rigidity parameter and Poisson’s ratio on the natural 
frequencies of an orthotropic rotating disc has been studied. 

The following conclusions are arrived at: 

1. As the clamping radius increases frequency increases. 

2. As the disc becomes thicker at the hub, frequency increases. 

3* As the disc becomes heavier towards periphery the 

frequency increases or decreases depending upon the 
rigidity parameter and disc configuration. 

4 . Increase in the rotational speed leads to the increase 
in frequency. 

5- If the ratio of Young's modulus in tangential 

direction to the Young's modulus in r°dial direction 
increases, frequency increases. 

6. Increases in Poisson's ratio lead to an increase or 
decrease in frequency depending on the value of the 
rigidity parameter. 
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Flexural natural f re^uencies^f o-r . the - t hr e-e -mod 
comparison with C.D. Mote, Jr. 


Radius 

a/b 

Nodal diameter 
n 

Frequency, w 

,f(pb 4 /SH 2 ) 

Present worm 

0 . D. Mote's wo rk 

0.1 

0 

2. 56 59 

0 

2. 5706 

0.2 

0 

3.1359 

3.7373 

0.3 

0 

4.0310 

4.0314 

0.4 

0 

5.4595 

5.4597 

0.5 

0 

7.8826 

7.8828 

0.6 

0 

12.4203 

12.4200 

0,7 

0 

22,3646 

22.3620 

0.8 

0 

51.1413 

50.9940 

0.1 

1 

2.1099 

2.1211 

0.2 

1 

2.9135 

2.9164 

0.3 

1 

3.9657 

3.9676 

0.4 

1 

5.5169 

5. 5181 

0.5 

1 

8.0433 

8 . 0447 

0.6 

1 

12.6710 

12.6720 

0.7 

1 

22.6949 

22.6950 

0.8 

1 

51.5420 

51.3980 


Continued 


Table 1 continued 


Iteidius 

a/b 

Nodal diameter 
n 

Fre^uencv, to 
- * n 

VC p b 4 /SH 2 ) 

Present work 

O.D. Mote’s work 

0.1 

2 

3.4077 

3.4147 

0.2 

2 

3.9036 

3.9196 

0.3 

2 

4.8158 

4.8367 

0.4 

2 

6.3335 

6.3452 

0.5 

2 

8.8991 

8.9051 

0.6 

2 

13.6184 

13.6220 

0.7 

2 

23.7714 

23.7750 

o 

• 

00 

2 

52.7713 

52.6 400 

oci 

3 

7.5361 

7.5338 

0.2 

3 

7.6360 

7.6509 

0.3 

3 

8.0355 

8.0644 

0.4 

3 

9.0562 • 

9.0806 

0.5 

3 

11.2342 

11.2500 

0.6 

3 

15.7238 

15.7330 

0.7 

3 

25.8156 

25.8250 

0.8 

3 

54.9052 

54.7940 


Note* H s i peripheral disc thickness. 

C. 

i.i »■ 'mmjr 

A 62179 
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Tabl e 2 

Comparison of approximate flexural natural frequencies with, 
exact flexural natural frequencies - Southwell. 


Radius 

a/b 

Nodal diameter 
n 

Frequency, w 

V( Pb 4 /3H 2 ) 

Present work 

Southwell's work 

0.276 

0 

3.778 

3.78 

0.642 

0 

15.581 

15.10 

0.840 

0 

80.49 

49.10 

0.060 

1 

1.808 

1.71 

0.597 

1 

5.459 

5.45 

0.603 

1 

12,866 

12.80 

0,634 

1 

15.172 

15.10 

0.771 

1 

39.195 

38.80 

0.827 

1 

69.100 

73.20 

0.186 

2 

3.812 

3.78 

0.394 

2 

6.219 

5.45 

0.522 

2 

9.690 

9.70 

0.769 

2 

39.692 

38.80 

0.810 

2 

58. 419 

60.70 

0.430 

3 

9.553 

9.71 

0.590 

3 

15.108 

15.15 

0.710 

3 

27.444 

29.70 

0.810 

3 

60.568 

60.60 

Notes 

H = Peripheral 

disc thickness. 
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gable 3 

Plexural natural frequencies - Variable thickness disc 
Ch(r) = H(r/b)-°- 9 ), Comparison with O.B. note, Jr. 


Badius 

a/b 

Nodal diameter 
n 

Frequency 

cfe V(P b 4 /M 2 ) 

Present work 

C.D.Mote’s work 

0.1 

0 

9.2363 

9.2499 

0.2 

0 

9.5081 

9.5135 

0.3 

0 

10.2617 

10.2650 

0.4 

0 

11.7753 

11.7770 

0.5 

0 

14.5742 

14.5760 

0.6 

0 

19.9915 

19.1330 

0.7 

0 

31.4183 

31.4160 

0.8 

0 

63.5736 

63.3700 

0.1 

1 

9.0925 

9.1027 

0.2 

1 

9.4598 

9.4631 

0.3 

1 

10.3140 

10.3160 

0.4 

1 

11.9149 

11.9170 

0.5 

1 

14.7880 

14.7900 

0.6 

1 

20.1912 

20.1940 

0,7 

1 

31.7594 

31.7590 

0,8 

1 

63.9740 

63.7650 


Continued. . . 
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labl e 3 c ooti'.'t inu ed . 


Radius 

aA 

WooSal diameter 
n 

Frequency, as 

V( p vVbj 2 ) 

Present work 

O.D.Mo'te’s wori 

0.1 

2 

9.7168 

9.7230 

0.2 

2 

10.1608 

10,1640 

0.3 

2 

11.0862 

11.0890 

0.4 

2 

12.7527 

12.7550 

0.5 

2 

15.6957 

15.6990 

0.6 

2 

21,1829 

21,1*90 

0.7 

2 

32.8564 

32 .8580 

0.8 

2 

65.2010 

6 4.9690 

— — - - - 




0.1 

3 

13.0749 

13 .1000 

0.2 

3 

13.2705 

13.2750 

0.3 

3 

13.8798 

13.8860 

0.4 

3 

15.2408 

15.2 470 

0. 5 

3 

17.9402 

17.9490 

0. 6 

3 

23.2756 

23.2890 

0.7 

3 

34.9070 

34.9120 

0. 8 

3 

67.3277 

67.0250 


F oto&JSi (i) E = 4r Peripheral disc thickness. 

, - (ii) All frequencies are for the 

zero nodal circle mode of vibration. « 
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lable 4 

Frequencies of Orfchotrmvi ,» - 

-opic su-tionary discs of constant 

thickness and constant density. 

Comparison with fiamaiah and Vi jay Sumar. 

(Zero nodal circle and 1 nodal diameter). 



Clamping radius 
a/b 

Value of 
constant 

M 

Frequency, 

F resent work 

( Vo “n 

# 

From appro::, formula" 
(.Bamaiah and Vijajr 
Kothp tO 

h25 

0.1 

0.548 

1.916 

.... Cv JL / 

1.931 

'.25 

0.3 

0.312 

3.530 

3.528 

.25 

0.5 

0.140 

7.461 

7.457 

.25 

0.7 

0.42 

22.000 

22.009 

.33 

0.1 

0.578 

1.998 

2.011 

.33 

0.3 

0. 312 

3.589 

3.587 

.33 

0.5 

0.140 

7.509 

7.505 

.33 

0.7 

O.O 42 

22.041 

22.049 

.50 

0.1 

0.548 

2.154 

2.163 

.50 

0.3 

0.312 

3.705 

3.703 

.50 

0.5 

0.140 

7.604 

7.601 

.50 

0.7 

0.042 

22.122 

22.128 

.00 

0.1 

0.578 

2.565 

2.565 

,00 

0.3 

0.312 

4.031 

4.031 

,00 

0.5 

0.140 

7.882 

7.882 

,00 

0.7 

O.O 42 

22.364 

22.364 



Frequency 


MOTE AND PRESENT 


I SOUTHWELL 



C lamping radius , a / b 

IG. 2 VARIATION OF FREQUENCY WITH CLAMPING RADIUS 

(COMPARISON WITH SOUTHWELL) 













STRESS DISTRIBUTION WITH RADIUS ( tfg- 0 3 , a I b = 0-5 

( 3 = 0 - 0 ). 

Darameter . 







Fig. 8 Variation of frequency with density parameter. 




10 VARIATION OF FREQUENCY WITH ROTATION SPEED 
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VARIATION OF FREQUENCY WITH ROTATION SPEED 








APPENDIX 


Approximate formula for determination of natural 
frequency of an orthotropic circular Plate : 

Ramaih and Vijay Kumar [ 22 '] has given an approximate 
formula. Natural frequency of an orthotropic plate can be 
determined if the frequency of the similar isotropic plate is 
known, from the following relation; 


0) 


n 


= f 1-M (1- 


D. 


) 3 


0) 


nx 


D 


d 9 < 


( A— 1 ) 


and 


or 


n 


D_ 


D 


e 


= [ i- m’ Cl — ^ ) 3 


e 


CO' 


nx 


D 




( A-2 ) 


where, n = circular frequency in radians per second 

ni = circular frequency in radians per second in 
isotropic case 

and the constants M and M are related by the equation 

M =1 (.1-3) 


The variation of M with clamping radius a/b is 


given in the following table. 



IX 



0.1 

0.578 

- 

0.413 

0.2 

0.436 

- 

0. 242 

0.3 

0.312 

- 

0.105 

0.4 

0.212 

0.100 

0. 036 

0.5 

0.140 

0.075 

0.019 

0.6 

0.080 

0. 050 

0.009 

0.7 

0.042 

0.029 

0. 005 

0.8 

0. 018 

0. 012 

0.002 
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